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ABSTRACT 

We discuss the eigenvalue problem for 2x2 and 3x3 octonionic Hermitian 
matrices. In both cases, we give the general solution for real eigenvalues, and 
we show there are also solutions with non-real eigenvalues. 



1. INTRODUCTION 

Finding the eigenvalues and eigenvectors of a given matrix is one of the basic techniques in 
linear algebra, with countless applications. The simplest case is that of (complex) Hermitian 
matrices, generalizing the familiar case of (real) symmetric matrices. This simple case is 
nevertheless very important, for instance in quantum mechanics, where the fact that such 
matrices have real eigenvalues allows them to represent physically observable quantities. 

The eigenvalue problem is usually formulated over a field, typically either the real num- 
bers R or the complex numbers C. We consider here the generalization to the other normed 
division algebras, namely the quaternions EI and the octonions O. We find that most of the 
basic properties are retained, provided they are reinterpreted to take into account the lack 
of commutativity of EI and O, and the lack of associativity of O. 

Our main result is the solution of the real eigenvalue problem for 3x3 octonionic 
Hermitian matrices, also known as Jordan matrices. It is straightforward to show [1] that 
such matrices admit 24 real eigenvalues, corresponding to eigenvectors which are independent 
over M. We show that these eigenvalues do not satisfy the characteristic equation even 
though the matrix itself does. Instead they generically come in 6 sets of multiplicity 4 
rather than the expected 3 sets of multiplicity 8 1 . We further show how to generalize the 
notion of orthogonality to the nonassociative case, recovering the standard decomposition of 
a Hermitian matrix in terms of its eigenvalues and eigenvectors. 

We begin in Section 2 with a review of the standard eigenvalue problem for real and 
complex Hermitian matrices, and then consider the quaternionic eigenvalue problem in Sec- 
tion 3. A brief discussion of the properties of octonions and octonionic matrices appears in 



1 After this work was completed, we discovered the existence of earlier work (in Russian) by Ogievetskh [2] 
which also claims 6 real eigenvalues for such matrices. 
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Section 4, after which the octonionic eigenvalue problem is considered for 2x2 and 3x3 
octonionic Hermitian matrices in Section 5 and Section 6, respectively. Finally, we discuss 
our results in Section 7. 

2. THE STANDARD EIGENVALUE PROBLEM 

The eigenvalue problem as usually stated is to find solutions A, v to the equation 

Av = Xv (1) 

for a given square matrix A. The basic properties of the eigenvalue problem for nxn complex 
Hermitian matrices are well-understood. 2 

Lemma 1^: An nxn complex Hermitian matrix A has n real eigenvalues 
(counting multiplicity). 

Proof: We give here only the proof that the eigenvalues are real. Let A, 
v, A satisfy (1), with = A. Then 

Xv^v = (Av)^v = v^Av = Xv^v (2) 

so that if v 7^ we have v^v ^ 0, which forces A = A. ■ 

Lemma 2^: Eigenvectors of an nxn complex Hermitian matrix A correspond- 
ing to different eigenvalues are orthogonal. 

Proof: For m — 1, 2, let v m be an eigenvector of A = A^ with eigenvalue 
X m . By the previous lemma, A m G K. Then 

\\v\v2 = (Avi)^V2 = v\Av2 = \2v\v2 (3) 

Then either Ai = A2 or v\v2 = 0. ■ 

Lemma 3^: For any nxn complex Hermitian matrix A, there exists an or- 
thonormal basis of C n consisting of eigenvectors of A. 

Proof: If all eigenvalues have multiplicity one, the result follows from the 
previous lemma. The Gram-Schmidt orthogonalization process 
can be used on any eigenspace corresponding to an eigenvalue 
with multiplicity greater than one. ■ 

These lemmas are equivalent to the standard result that a complex Hermitian matrix can 
always be diagonalized by a unitary transformation. It is important for what follows to 
realize that the form of the proofs given above relies on both the commutativity and the 
associativity of C. 

Combining the above results, it is easy to see that any (complex) Hermitian matrix A 
admits a decomposition in terms of an orthonormal basis of eigenvectors. 



2 We could just as well have started with the case of real symmetric matrices. 
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Theorem 1^: Let A be an n x n complex Hermitian matrix. Then A can be 
expanded as 

n 

A = ^mVmvtn (4) 
m=l 

where {v m ; m — 1, ...,n} is an orthonormal basis of eigenvectors 
corresponding to eigenvalues X r 



Proof: By the previous lemma, there exists an orthonormal basis {v m } 
of eigenvectors. It then suffices to check that 

n 

J2 ^mVmvlnVk = \ k V k (5) 

771=1 

But this follows by direct computation using ort honor mality. ■ 

Furthermore, the set of eigenvalues {A m } is unique, and the (unit) eigenvectors are unique 
up to unitary transformations in the separate eigenspaces (which reduce to multiplication 
by a complex phase for eigenvalues of multiplicity one). 

3. THE QUATERNIONIC EIGENVALUE PROBLEM 

The quaternions EI double the dimension of the complex numbers by adding two addi- 
tional square roots of —1, usually denoted j and k. The multiplication table follows from 

i 2 = j 2 = k 2 = —1 ij = k = —ji 

and associativity; note that EI is not commutative. Equivalently, H can be viewed via the 
Cayley-Dickson process as the sum of 2 copies of the complex numbers 

H = C + Cj (6) 

with k being defined by k = ij. 

The eigenvalue problem (1) for Hermitian matrices A over EI immediately yields the first 
unexpected result: The eigenvalues need not be real. An example is given by 

Furthermore, because of the lack of commutativity, multiples of eigenvectors are not neces- 
sarily eigenvectors. For instance, the vector 



is an eigenvector of the matrix 



* = li-°, 1+ o'') 



with eigenvalue v^2, but jv\ is not an eigenvector of A\. This example illustrates an impor- 
tant point: We must distinguish between right and left multiplication. Since 

A(vq) = (Av)q (10) 
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by associativity, right multiples of eigenvectors are indeed eigenvectors. For example, v\j is 
an eigenvector of the matrix A\ above, with the same eigenvalue (y/2). 

Similarly, we must carefully distinguish between the left eigenvalue problem (1) and the 
right eigenvalue problem 

Av = v\ (11) 
It turns out that all the right eigenvalues of (quaternionic) Hermitian matrices are real. 

Lemma ljj: The right eigenvalues of annxn quaternionic Hermitian matrix 
A — A' are real. 

Proof: This is just a careful rewrite of (2). Explicitly, 

\(v*v) = (\ v i) v = (Av)^v = (v^A)v = v\Av) = v^(v\) = (y*v)\ (12) 
and the result follows since v*v e K. ■ 

In the above, A is now the quaternionic conjugate of A and f denotes (quaternionic) Hermitian 
conjugation. Once the eigenvalues have been shown to be real, orthogonality of eigenvectors 
with different eigenvalues follows as in the complex case. 

Lemma 2jj: Right eigenvectors of an n x n quaternionic Hermitian matrix A 
corresponding to different eigenvalues are orthogonal. 

Proof: This is just a rewrite of (3). Explicitly, 

X\{v\v 2 ) = {Xiv\)v 2 = (Av 1 )^v 2 = {v\A)v 2 

t t t ^ ' 

= v[(Av 2 ) = v[(v 2 \ 2 ) = (v[v 2 )X 2 

so that, since X m G K, either Ai = A2 or v\v 2 = 0. ■ 

Finally, since the (right) eigenvalues are real, (right) multiples of eigenvectors are still eigen- 
vectors. 

Putting it all together, we obtain a decomposition of any (quaternionic) Hermitian matrix 
A of the form (4), where the real eigenvalues {A m } and their eigenspaces are again unique. 

Lemma 3jj: For any n x n complex Hermitian matrix A, there exists an or- 
thonormal basis ofW 1 consisting of eigenvectors of A. 

Theorem ljj: Let A be annxn quaternionic Hermitian matrix. Then A can be 
expanded as in (4), where {v m ; m — 1, ...,n} is an orthonormal 
basis of eigenvectors of A, with real eigenvalues A m . 

The proofs of each of these results is identical to the complex case previously considered. For 
eigenvalues of multiplicity one, the (unit) eigenvectors are unique up to a (right) quaternionic 
phase. 

The right eigenvalue problem over EI is therefore just a straightforward extension of the 
complex eigenvalue problem [3,4,5]. The left eigenvalue problem turns out to be of consid- 
erable interest as well, and will be considered elsewhere in the context of 2 x 2 octonionic 
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Figure 1: The representation of the octonionic multiplication table using the 7-point projective 
plane. Each of the 7 oriented lines represents a quaternionic triple. 

Hermitian matrices [6] 3 . 

4. OCTONIONS AND OCTONIONIC MATRICES 
a) The Octonions 

The octonions O can be viewed via the Cayley-Dickson process as the direct sum of two 
copies of the quaternions [7] 

= M + M£ (14) 

where £ is yet another square root of —1. The octonions are thus spanned by the identity 
element 1 and the 7 imaginary units {i,j,k,£,i£,j£,k£}. These units can be grouped into 
associative "triples" in 7 different ways, each of which generates (the imaginary part of) a 
quaternionic subspace. Any three of these imaginary units which do not lie in a such a triple 
anti-associate. The multiplication table can be neatly summarized by appropriately labeling 
the 7-point projective plane, as shown in Figure 1. For a good introduction to the octonions, 
including some applications to physics, see [8,9]. 

When working with small numbers of octonions, it is important to realize that simpli- 
fications take place by virtue of the automorphism group G2 of O. For instance, a single 

3 Another argument that the right eigenvalue problem is the natural one is based on their use in diag- 
onalizing a matrix. Cohn [5] considers quaternionic matrices A which are diagonalizable, in the sense 
that there exists an invertible matrix U and a diagonal matrix D such that 

A = UDIJ- 1 

But this is equivalent to 

AU = UD 

so that the columns of U are eigenvectors of A with right eigenvalues taken from the elements of D. We 
will return to this issue below. (Cohn uses the term left eigenvalue to describe the eigenvector problem 
for row vectors multiplied on the right by A; this is completely different from our use of the same term.) 
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octonion may be assumed without loss of generality to be complex, i.e. to lie in the complex 
subspace of O spanned by Similarly, a second octonion can be assumed to lie in the 

subspace spanned by {l,i,j}, and a third can be assumed to lie in the subspace spanned 
by j, k, £}. Only when four or more octonions are involved is it necessary to consider 
"generic" octonions, involving all the basis directions. 

The octonions are not associative. Nevertheless, since any 2 octonions lie in a quater- 
nionic subspace, products involving only 2 different octonions (and their octonionic conju- 
gates) do associate. For example, 

p(pq) = p 2 q (15) 
which is a weak form of associativity known as alternativity. 
The squared norm of an octonion a is given by 

\a\ 2 := aa (16) 
where a denotes the (octonionic) conjugate of a. The commutator of a and b is given by 

[a,b]:=ab-ba (17) 

the associator of a, b, c is given by 

[a,b,c] := (ab)c- a(bc) (18) 

and we use to denote the (octonionic) Hermitian conjugate of the matrix A. Both the 
commutator and the associator are purely imaginary, totally antisymmetric, and change sign 
if any one of their arguments is replaced by its conjugate. Another octonionic product with 
the latter two properties is given by 4 

<S>(a,b,c) = ^Re([a,b]c) (19) 

b) Octonionic matrices 

The lack of associativity complicates the treatment of matrices with octonionic entries. 
While matrix multiplication can be defined for matrices of arbitrary size in the usual way, 
only the 2x2 and 3x3 octonionic Hermitian matrices form Jordan algebras [12]. We 
therefore limit ourselves to these two cases. 

Any complex number x = x\ + ix2 can be written as the real matrix a;^ 2 )' 

which allows (complex) matrix multiplication to be rewritten in terms of real matrices. This 
process can be generalized to the quaternions, but fails for the octonions — as it must, 
since octonionic multiplication is not associative. But the complex number x can also be 

written as the real vector ) ' l eao ^ m S to a representation of (complex) multiplication as 

the product of a (real) matrix with a (real) vector, and this latter process does generalize to 
the octonions. 



4 $(a, b, c) is in fact the same as the associative 3-form [10,11] 

$(a, b, c) = Rc(a x b x c) = ^Rc (a(bc) - c(bafj 

which reduces to the vector triple product when a, b, c are imaginary quaternions. 

(Note that a x b x c := \ (a(bc) — c(ba)) is the triple cross product, not the iterated cross product.) 
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A similar process can be used to represent matrices over one of these division algebra 
as matrices over any smaller such algebra. Under this transformation a complex Hermitian 
matrix is mapped to a symmetric real matrix, and quaternionic and octonionic Hermitian 
matrices can be transformed into either real symmetric matrices or complex Hermitian ma- 
trices. For example, an n x n octonionic Hermitian matrix can be mapped to an 8n x 8n 
symmetric real matrix (or a 4n x 4n complex Hermitian matrix). It would seem as if we 
could therefore reduce the eigenvalue problem to the real case, but this is misleading for 
several reasons. 

First of all, while the real formalism guarantees the existence of real eigenvalues, it does 
not rule out the possibility that eigenvalues might exist which are not real. This is because 
the general eigenvalue problem over a normed division algebra transforms into a matrix 
equation, which only reduces to an ordinary eigenvalue problem for real eigenvalues. As we 
will see, octonionic Hermitian matrices admit (right) eigenvalues which are not real. 

Secondly, in the real formalism it is not very easy to determine the multiplicity of the 
real eigenvalues. One might expect n octonionic eigenvectors with at most n different real 
eigenvalues. The real formalism does guarantee us 8n independent (over R) eigenvectors with 
real eigenvalues, but we are not guaranteed that the 8n eigenvalues occur with multiplicity 
8 (or a multiple thereof). In fact, we will see below that this is not the case. 

The final drawback of this approach is that the orthogonality of eigenvectors with differ- 
ent eigenvalues does not follow from the corresponding statement on the real, transformed 
eigenvectors. This is because the transformation does not preserve the inner product, which 
is real in one case and not in the other. Nevertheless, the real parts of the inner products 
do agree, so that 

v ] w + w ] v = (20) 

As we will see below, there is a stronger orthogonality condition on eigenvectors with different 
eigenvalues, which generalizes the usual notion of orthogonality. 

For all of these reasons, we choose to work directly with the octonionic matrices. 

5. 2 x 2 OCTONIONIC HERMITIAN MATRICES 

Are the eigenvalues of octonionic Hermitian matrices real? Consider first the special case 
where A and v (but not necessarily A) in (11) lie in a quaternionic subspace. Then (12) still 
holds, and alternativity allows us to shift the parentheses and conclude that A is real. 

For instance, since the diagonal entries of a 2 x 2 Hermitian matrix A are real, its com- 
ponents always lie in a complex subspace of O, and therefore A clearly possesses eigenvectors 
which lie in the same complex subspace, and which have real eigenvalues. Multiplication 
(on the right) of these eigenvectors by an arbitrary octonion q leads to a new eigenvector, 
which lies in the quaternionic subspace spanned by the single octonionic direction in A and 
the octonionic multiple q. Furthermore, this new eigenvector has the same real eigenvalue 
as the original eigenvector, that is 

Av = vX A{vq) = (Av)q = (v\)q = (vq)X (21) 

since A and v are complex and A is real. 

In general, however, the key use of associativity in the middle of the derivation of (12) 
is not allowed. We are thus led to suspect that there exist octonionic Hermitian matrices 
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which admit (right) eigenvalues which are not real. This turns out to be correct, as is shown 
by the following example: 



1 i\(j\_(j + it 
-i 1 ) U I U - k 



[1 - k£) (22) 



Further details for the case of octonionic Hermitian matrices whose eigenvalues are not 
real will appear elsewhere [6]. 

a) The real eigenvalue problem 

We now turn to the case of real eigenvalues. The general 2x2 octonionic Hermitian 
matrix can be written 

A=(P a ) (23) 

with p,m el and a G O, and satisfies its characteristic equation 

A 2 - (tr A)A+ (det A) I = (24) 

where tr A denotes the trace of A, and where there is no difficulty defining the determinant 
of A as usual via 

det A = pm - \a\ 2 (25) 
If a = the eigenvalue problem is trivial, so we assume a ^ 0. If we set 



then (11) leads to 



o={y) < 26 > 

from which it follows (unless v — 0) that 

det(A7 - A) = A 2 - (tr A) X + det A = (28) 
as usual. Eigenvectors can thus be written in either of the forms 



. a(A — p) 

where £ G O is arbitrary and where A is either of the 2 solutions of (28). 5 This shows 
that all eigenvectors of 2 x 2 Hermitian matrices with real eigenvalues are obtained from the 
usual complex eigenvectors by (right) multiplication by an arbitrary octonion. 

Lemma 1<q : Let A be a 2 x 2 complex Hermitian matrix. Then w is an oc- 
tonionic eigenvector of A with real eigenvalue A if and only if 
w = i>£ ; where £ G O is arbitrary and where v is a complex 
eigenvector of A with the same eigenvalue. 

Since a 2 x 2 octonionic Hermitian matrix contains only one independent octonion, we can 
assume without loss generality that any such matrix is complex, and thus apply this lemma 
to it. 



5 Note that x ^ ^ y if a ^ 0, and that in this case the 2 solutions of (28) are distinct. 
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b) Orthogonality and Decompositions 

As already noted, since A lies in a complex subspace of O, it admits a complete set 
of complex eigenvectors with real eigenvalues, which can be used to obtain the decompo- 
sition (4). But what about a decomposition in terms of the general solution of the real 
eigenvalue problem? 

The general solution is given by (29). Choosing the first form, we obtain a complete set 
of eigenvectors by considering both solutions A± to (28), obtaining 

« 3o) 

Are these eigenvectors orthogonal? Direct computation yields 

(31) 

where we have used 

a[a, x,y] = —a[a,y,x] = [a,ay,x] = [x,a,ay] = (xa)(ay) — \a\ 2 xy (32) 
or equivalently 

a((ax)y) = (xa)(ay) (33) 

for any octonions a, x, y. Thus, the eigenvectors are not necessarily orthogonal in the 
traditional sense except for the quaternionic eigenvalue problem, when the associator auto- 
matically vanishes. 

At first sight, this apparent lack of orthogonality seems to rule out a decomposition of 
the form (4). However, due to the lack of associativity, what is needed for (4) to hold is an 
appropriately generalized notion of orthogonality, namely 

Definition: Let v and w be two octonionic vectors. We will say that w is 
orthogonal to v if 

( OT t) w = o (34) 

The vectors {v, w} are orthonormal if in addition v'v = 1 = w^w. 

Direct computation shows that the eigenvectors v± above are indeed mutually orthogonal 
in this sense, which provides a computational proof of the following lemma. 

Lemma 2q 2 : If v and w are eigenvectors of the 2x2 octonionic Hermitian 
matrix A corresponding to different real eigenvalues, then v and 
w are mutually orthogonal in the sense of (34). 

Proof: From Lemma 1@ 2 , we can write 

v = va, w = v(3 (v, w G C; a, (3 G O) 

where C C O is the complex subspace containing the elements of 
A. But then vv* = \a\ 2 vv^ , and 

(vv^)w = \a\ 2 (vv^)(w/3) 

which associates since only 2 octonionic directions are involved. 
But v^w = by the usual properties of complex eigenvectors. ■ 



v' + V- 



— a a 
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In order for a decomposition of the form (4) to exist, we also need a vector version of 
alternativity, which in fact holds for octonionic vectors of any size: 6 

Proposition 1: For any octonionic vector v e O n , 

(vv^) v = v (v^v) (35) 

This proposition shows in particular that any normalized vector v is an eigenvector of the 
matrix vv* with eigenvalue 1, as required by (4). We conclude that the decomposition (4) 
holds unchanged for real eigenvalues. We thus have: 

Theorem 1<q : Let A be a 2 x 2 octonionic Hermitian matrix. Then A can be 
expanded as in (4), where {fi,t>2} are orthonormal (as per (34) ) 
eigenvectors of A corresponding to the real eigenvalues X m . 

Proof: Provided the real eigenvalues of A are distinct, Lemma 2q guar- 
antees the existence of orthonormal eigenvectors, which are also 
eigenvectors of the decomposition (4) with the same eigenvalues, 
and the result follows. But if A has a repeated eigenvalue, it 
must be a multiple of the identity matrix, for which the result 
holds. ■ 

Using the same technique as in Lemma 2q , it is straightforward to show that 

(vv*)(vv*) = (v*v)(vv') (36) 

for any v E O 2 , and that 

(vv^)(ww f ) = (37) 

if v and w are eigenvectors of A with distinct real eigenvalues, since each term in parentheses 
lies in C. The decomposition in the preceding theorem is thus in terms of orthonormal 
idempotents V{v\. We also have 

^ + ^ = 7 (38) 

which could be proved directly using the fact that the left-hand side has repeated eigen- 
value 1. Furthermore, since by Lemma 2q 2 A and v contain only 2 octonionic directions, 

(Av)v^ = A(vv ] ) (39) 
which leads to the following one-line alternative derivation of Theorem 1q 

A = A[j2 Viv}] = E A( Vi vj) = j2(A Vi )vj = E Wj (40) 

\i=l ) i=l i=l i=l 

As we will see below, however, this argument relies heavily on Lemma 2q 2 , which fails in 
the 3x3 case. 

We have not yet discussed whether the orthonormal eigenvectors in the preceding theo- 
rem constitute a basis of O 2 . For any orthonormal vectors, we have the following lemma. 



For 2-component vectors, this proposition is just the 3-\&'s rule of supersymmetry theory [13]. 
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Lemma 3q 2 : Let v,w E O 2 be orthonormal in the sense of (34), and let g be 
any vector in O 2 . Then 

g = (W) g + ( ww t) g 

Proof: This follows immediately from (38). ■ 

In the associative case, this lemma shows how to write any vector g in terms of its components 
along v and w, thus establishing {v, w} as a basis. One could adopt similar language in the 
nonassociative case, although the "component" of g "along" v would no longer point in the 
v direction. Nevertheless, this terminology is extremely attractive, as it allows the Gram- 
Schmidt orthogonalization process to be used to determine the component of one vector 
orthogonal to another. 

Proposition 2: Let v,w £ O 2 . Then 

(W) ( w _ (™1 w \ = o (41) 

\ V'V J 

Proof: This follows from the alternativity of 2 x 2 octonionic Hermitian 
matrices and (36). ■ 

6. 3x3 OCTONIONIC HERMITIAN MATRICES 

We now turn to the 3x3 case. It is not immediately obvious that 3x3 octonionic 
Hermitian matrices have a well-defined determinant, let alone a characteristic equation. 
We therefore first review some of the properties of these matrices before turning to the 
eigenvalue problem. As in the 2x2 case, over the octonions there will be solutions of 
the eigenvalue problem with eigenvalues which are not real; we consider here only the real 
eigenvalue problem. 

a) Jordan matrices 

The 3x3 octonionic Hermitian matrices, henceforth referred to as Jordan matrices, form 
the exceptional Jordan algebra (also called the Albert algebra) under the Jordan product 7 

AoB := ^(AB + BA) (42) 

which is commutative, but not associative. A special case of this is 

A 2 = A o A (43) 

and we define 

A 3 := A 2 o A = A o A 2 (44) 

Remarkably, with these definitions, Jordan matrices satisfy the usual characteristic equa- 
tion [10] 

A 3 - (tr A) A 2 + a(A) A - (det A) I = (45) 

where cr(A) is defined by 

a(A) := i((tr A) 2 - tr (A 2 )) (46) 



7 The 2x2 octonionic Hermitian matrices form a special Jordan algebra since they are alternative [12]. 
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and where the determinant of A is defined abstractly in terms of the Freudenthal product. 8 
Concretely, if 

(p a b\ 
a m c (47) 
b c n J 

with p,m,n G M and a, b, c e O then 

tr A = p + m + n 

a (A) =pm + pn + mn-\a\ 2 - \b\ 2 - |c| 2 (48) 



det A = pmn + b(ac) + b(ac) — n\a\ 2 — m\b\ z — p\c\' z 

b) The real eigenvalue problem 

As discussed above, n x n Hermitian matrices over any of the normed division algebras 
can be rewritten as symmetric kn x kn real matrices, where k denotes the dimension of the 
underlying division algebra, it is clear that a 3 x 3 octonionic Hermitian matrix must have 
8 x 3 = 24 real eigenvalues [1]. However, as we now show, instead of having (a maximum 
of) 3 distinct real eigenvalues, each with multiplicity 8, we show that there are (a maximum 
of) 6 distinct real eigenvalues, each with multiplicity 4. 

The reason for this is that, somewhat surprisingly, a (real) eigenvalue A of a Jordan 
matrix A does not in general satisfy the characteristic equation (45). 9 To see this, consider 
the eigenvalue equation (1), with A as in (47), AgI, and where 

(49) 

\z J 

Explicit computation yields 

(A — p)x = ay + bz (50) 
(X — m)y = cz + ax (51) 
(A - n)z = bx + cy (52) 

so that 

(A — p)(X — m)y = (A — p)(cz + ax) = (A — p)cz + a(ay + bz) (53) 



8 The Freudenthal product of two Jordan matrices A and B is given by [14] 

A * B = A o B - i (A tr (B) + B tr (A)j + ^ (tr (A) tr (B) - tr (A o B)j 



The determinant can then be defined as 



det (A) = itr ({A * A) o A 



9 Ogievetskii [2] constructed a 6th order polynomial satisfied by the real eigenvalues, which he called 
the characteristic equation. This polynomial is presumably equivalent to the modified characteristic 
equations (for both values of r) given below. 
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which implies 

(A — p)(X — m) — \a\ 2 y — a(bz) + (A — p)cz 
Assume first that A ^ p. Using (50) and (54) in (52) leads to 



(54) 



(X-p)(X-m) 



(A — p)(X — n)z 



(A — p)(X — m) 





2 


\a 






2 


\a 





(A — p)(bx + cy) 



(A — p)(X — m) — \a\ b(ay + bz) + (A — p)c a(bz) + (A — p)cz 



= b 



a (a(bz) + (A — p)cz} + (A— p)(X — m) 



b{bz) 



(55) 



(A-P) 



+ (A — p)c a(bz) + (A — p)cz 
(X-m)\b\ 2 z + (A - p)\c\ 2 : z + b(a(czf) +c(a(h 



Expanding this out and comparing with (48) results finally in 
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[ det (XI-A)]z= [A 3 - (trA) A 2 + a(A) X - det A 



b[a(cz) 



c\a[bz 



b(ac) + (ca)b 



Now consider the case A = p. We still have (54), which here takes the form 

— \a\ 2 y = a(bz) 

Inserting this into (51), we can solve for x, obtaining 

— |a| 2 x = a(cz) + (p — m)bz 
Finally, inserting (57) and (58) in (52) yields 



(Ja| 2 (p — n) + \b\ 2 (p — mfj z = b(a(cz)^j + c (a(bz)^) 



(56) 

(57) 
(58) 
(59) 



(60) 



Comparing with (48) and using A = p, we see that (56) still holds, and thus holds in general. 

If a, b, c, and z associate, the RHS of (56) vanishes, and A does indeed satisfy the 
characteristic equation (45); this will not happen in general. However, since the LHS of (56) 
is a real multiple of z, this must also be true of the RHS, so that 

b(a(cz)^J + c (a(bzfj — b(ac) + (ca)b z = rz r£l 

which can be solved to yield a quadratic equation for r as well as constraints on z. 

Lemma lo 3 : The real eigenvalues of the 3x3 octonionic Hermitian matrix A 
satisfy the modified characteristic equation 

det(A/ - A) = A 3 - (tr A) X 2 + a{A) X - det A = r 

where r is either of the two roots of 



r 2 + 4$(a, b,c)r - [a, 6, c] 







(61) 
(62) 



with a, b, c as defined by (47) and where $ was defined in (19). 
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We have recently been able to generalize this to the case where A is not real [6]. 

13 



Proof: These results were obtained using Mathematica to solve (60) by 
brute force for real r and octonionic z given generic octonions 
a, b, c [15]. ■ 



Furthermore, provided that [a, b, c] 7^ 0, each of x, y, and z can be shown to admit an 
expansion in terms of 4 real parameters. 

Corollary 1: With A and r as above, and assuming [a, b, c] ^ 0, 



The real paramaters a, (3, 7, 5 may be freely specified for one (nonzero) component, say z; the 
remaining components x,y have a similar form which is then fully determined by (50)-(52). 

Corollary 2: The real eigenvalues of A are the same as those of A. 

Proof: Direct computation (or (48)) shows that 



The solutions of (61) are real, since the corresponding 24 x 24 real symmetric matrix has 
24 real eigenvalues. We will refer to the 3 real solutions of (61) corresponding to a single 
value of r as a family of eigenvalues of A. There are thus 2 families of real eigenvalues, each 
corresponding to 4 independent (over M) eigenvectors. 

We note several intriguing properties of these results. If A is in fact complex, then the 
only solution of (62) is r = 0, and we recover the usual characteristic equation with a unique 
set of 3 (real) eigenvalues. If A is quaternionic, then one solution of (62) is r = 0, leading 
to the standard set of 3 real eigenvalues and their corresponding quaternionic eigenvectors. 
However, unless a, b, c involve only two independent imaginary quaternionic directions (in 
which case $(a, b,c) — — [a, b,c]), there will also be a nonzero solution for r, leading to a 
second set of 3 real eigenvalues. From the preceding corollary, we see that this second set 
of eigenvalues consists precisely of the usual (r = 0) eigenvalues of A\ Furthermore, since 
A(£v) = £(Av) for A,ti £ I and imaginary f GO orthogonal to H, the eigenvectors of A 
corresponding to r 7^ are precisely £ times the quaternionic (r = 0) eigenvectors of A. In this 
sense, the octonionic eigenvalue problem for quaternionic A is equivalent to the quaternionic 
eigenvalue problem for both A and A together. Finally, if A is octonionic (so that in particular 
[a, b, c] 7^ 0), then there are two distinct solutions for r, and hence two different sets of real 
eigenvalues, with corresponding eigenvectors. Note that if det A = 7^ [a, b, c] then all of 
the eigenvalues of A will be nonzero! 

The final suprise lies with the orthogonality condition for eigenvectors v, w corresponding 
to different eigenvalues. It is not true (cf. (31)) that v^w = 0, although the real part (20) of 




(63) 



with a, /3, 7, 5 e M. Similar expansions hold for x and y. 



det A = det A - 4$(a, b,c) 



(64) 



But — 4$(a, b, c) is precisely the sum of the roots of (62), and 
replacing A by A merely flips the sign of r, that is r[A] = —r[A]. 
Thus, the 2 possible values of det A + r[A] are precisely the same 
as those for det A + r[A]. Since tiA = tr A and a (A) = a (A), 
(61) is unchanged. ■ 
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this expression does vanish. However, just as in the 2x2 case, what is needed to ensure a 
decomposition of the form (4) is (34), and a lengthy, direct computation verifies that (34) 
holds provided that both eigenvectors correspond to the same value of r. 

Lemma 2q 3 : If v and w are eigenvectors of the 3x3 octonionic Hermitian 
matrix A corresponding to different real eigenvalues in the same 
family (same r value), then v and w are mutually orthogonal in 
the sense of (34). 

Proof: The modified characteristic equation (61) can be used to elimi- 
nate cubic and higher powers of A from any expression. Further- 
more, given two distinct eigenvalues Ai ^ A 2 , subtracting the two 
versions of (61) and factoring the result leads to the equation 

(A? + Ai A 2 + A|) - tr A(X 1 + A 2 ) + a(A) = (65) 

which can be used to eliminate quadratic terms in one of the 
eigenvalues. 11 ■ 

For Jordan matrices, we thus obtain two decompositions of the form (4), corresponding 
to the two sets of real eigenvalues. For each, the eigenvectors are fixed up to orthogonal 
transformations which preserve the form (63) of z. 

Theorem 1q 3 : Let A be a 3 x 3 octonionic Hermitian matrix. Then A can be 
expanded as in (4), where {vi,V2,vs} are orthonormal (as per 
(34) ) eigenvectors of A corresponding to the real eigenvalues \ m , 
which belong to the same family (same r value). 

Proof: Fix a family of real eigenvalues of A by fixing r. If the eigenvalues 
are distinct, then the previous theorem guarantees the existence 
of orthonormal eigenvectors, which are also eigenvectors of the 
decomposition (4) with the same eigenvalues, and the result fol- 
lows. 

If the eigenvalues are the same, the family consists of a sin- 
gle real eigenvalue A with multiplicity 3. Then tr (A) = 3A and 
cr (A) = 3A 2 . Writing out these two equations in terms of the 
components (47) of A, and inserting the first into the second, 
results in a quadratic equation for A; the discriminant D of this 
equation satisfies D < 0. But A is assumed to be real, which 
forces D — 0, which in turn forces A to be a multiple of the 
identity matrix, for which the result holds. 

The remaining case is when one eigenvalue, say //, has multi- 
plicity 2 and one has multiplicity 1. Letting v be a (normalized) 
eigenvector with eigenvalue fj,, consider the matrix 

X = A - a vv ] (66) 



We used Mathematica to implement these simplifications in a brute force verification of (34) in this 
context, which ran for 6 hours on a SUN Sparc20 with 224 Mb of RAM [15]. 
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with a G R. For most values of a, X will have 3 distinct real 
eigenvalues, whose eigenvectors will be orthogonal by the pre- 
vious theorem. But this means that eigenvectors of X are also 
eigenvectors of A; the required decomposition of A is obtained 
from that of X simply by solving (66) for A. ■ 

Note in particular that for some quaternionic matrices with determinant equal to zero, one 
and only one of these two decompositions will contain the eigenvalue zero. 

In the 2x2 case, (36) tells us that, for normalized v, vv** squares to itself, and hence 
is idempotent. As already noted, the decomposition of Theorem 1q is thus an idempotent 
decomposisition. But (36) fails in the 3x3 case, so that the decomposition in Theorem 1q 
is therefore not an idempotent decomposition. 

It is nevertheless straightforward to show, in analogy with (38), that if u, v, and w are 
orthonormal in the sense of (34), then 

uu* + vv* + wvo* = I (67) 

since the left-hand side has eigenvalue 1 with multiplicity 3. This permits us to view {u, v, w} 
as a basis of O 3 in the following sense 

Lemma 3q 3 : Let u,v,w G O 3 be orthonormal in the sense of (34), and let g 
be any vector in O 3 . Then 

g = {uu*) g + {vv*) g + {ww"*) g (68) 

Proof: This follows immediately from (67). ■ 

However, another consequence of the failure of (36) in the 3x3 case is that the Gram- 
Schmidt orthogonalization procedure (41) no longer works. It appears to be fortuitous that 
we are nevertheless able to find orthonormal eigenvectors in the 3x3 case with repeated 
eigenvalues; we suspect that this might fail in general, perhaps already in the 4x4 case with 
an eigenvalue of multiplicity 3. 

7. DISCUSSION 

Our main result is (61) together with (62), which shows how to modify the characteristic 
equation for a 3 x 3 octonionic Hermitian matrix in order to find its real eigenvalues, and 
which further shows that there are in general 2 families of solutions of these equations, each 
consisting of 3 eigenvalues of multiplicity 4. 

We have further shown how to use the corresponding eigenvectors to recover the usual 
decomposition (4) of a Hermitian matrix in terms of its eigenvectors for both 2x2 and 3x3 
octonionic Hermitian matrices. In the process we were led to introduce an appropriately 
generalized notion of orthogonality, namely (34). 

We can relate our notion of orthonormality to the usual one by noting that n vectors in 
O n which are orthonormal in the sense (34) satisfy 

vv* + ... + ww^ = I (69) 

If we define a matrix U whose columns are just v , w, then this statement is equivalent to 

Urf = I (70) 
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Over the quaternions, left matrix inverses are the same as right matrix inverses, and we 
would also have 

U ] U = I (71) 

or equivalently 

v*v = 1 = ... = w^w; v^w = = ... (72) 

which is just the standard notion of orthogonality. These two notions of orthogonality fail to 
be equivalent over the octonions; we have been led to view the former as more fundamental. 

We can now rewrite the eigenvalue equation (11) in the form 

AU = UD (73) 

where D is a diagonal matrix whose entries are the real eigenvalues. Multiplying (73) on the 
left by C/t yields 

U\AU) = U\UD) = (U ] U)D (74) 

(since D is real), but this does not lead to a diagonalization of A since, as noted above, 
is not in general equal to the identity matrix. However, Theorem 1q can be rewritten as 

A = UDU ] (75) 

so that in this sense A is diagonalizable. Furthermore, multiplication of (73) on the right by 
C/t shows that 

(AU)U ] = {UD)U ] = A = A{UU ] ) (76) 

and this assertion of associativity can be taken as a restatement of both Theorem 1q and 
Theorem 1q . In the 2x2 case, this associativity holds for a single eigenvector v, which is 
(39), and which was used in the one-line proof (40) of Theorem 1q . However, (39) fails in 
the 3x3 case, and we are unaware of a correspondingly elegant proof of Theorem 1q . 

Many of our results were obtained using Mathematica [15]. In fact, the only proof we 
currently have of the 3x3 orthogonality result, namely Lemma 2q 3 , uses Mathematica to 
explicitly perform a horrendous, but exact, algebraic computation. While one could hope for 
a more elegant mathematical proof of this result, the Mathematica computation nevertheless 
establishes a result which would otherwise remain for the moment merely a conjecture. This 
is a good example of being able to use the computer to verify one's intuition when it may 
not be possible to do so otherwise. This issue is further discussed in [15]. 

Finally, it is intriguing that (some) Hermitian octonionic matrices admit eigenvalues 
which are not real. In particular, this means that octonionic self-adjoint operators do not 
necessarily have a (purely) real spectrum. We plan to report separately on these matters [6]. 
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